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Abstract 

To analyze linear field equations on a locally homogeneous spacetime 
by means of separation of variables, it is necessary to set up appropriate 
harmonics according to its symmetry group. In this paper, the harmonics 
are presented for a spatially compactified Bianchi II cosmological model 
— the mlgeometnc model. Based on the group structure of the Bianchi II 
group (also known as the Heisenberg group) and the compactified spatial 
topology, the irreducible differential regular representations and the mul- 
tiplicity of each irreducible representation, as well as the explicit form of 
the harmonics are all completely determined. They are also extended to 
vector harmonics. It is demonstrated that the Klein-Gordon and Maxwell 
equations actually reduce to systems of ODEs, with an asymptotic solu- 
tion for a special case. 

1 Introduction 

A basic strategy to analyze linear field equations on a given spacetime, like 
linear perturbation equations of Einstein's equation, is to separate the equa- 
tions using appropriate harmonics. The harmonics for a given manifold are in 
general determined by the underlying symmetry group (isometry group) and 
the topology of the manifold. The simplest example is that of a commutative 
group acting on a flat manifold. If the manifold is compactified to, e.g., a torus, 
then functions on the manifold are expanded in the form of usual Fourier se- 
ries. When the group is noncommutative, however, the harmonics become much 
more complicated. The most familiar example in homogeneous cosmology is the 
SU(2) (Bianchi IX) case [21151^], where one needs to use the spherical har- 
monics (and their generalization to vector and tensor harmonics if necessary) to 
separate field variables. A notable nontrivial example is the x R (Bianchi 
III) case |2()l IT?i) with compactified three dimensional manifold. Since such a 
manifold is a direct product of two submanifolds, a closed hyperbolic plane and 
a circle, the harmonics are simply given by making products of those for the 
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two lower dimensional manifolds. Note that the Bianchi III belongs to Class 
B 1^. Separations of variables regarding locally rotationally symmetric (LRS) 
Class A Bianchi types (and gravitational perturbations) were discussed in 
without compactification. 

In this paper we consider the generic (i.e., non-LRS) Bianchi II model with 
compactification. The Bianchi II type is one of the class A types. The underly- 
ing symmetry group is the 3-dimensional Heisenberg group Hi , which we refer 
to as the Bianchi II group Gn. The Bianchi II homogeneous manifolds corre- 
spond to Thurston's nilgeometry |2S|- We will refer to spatially compactified 
Bianchi II spacetimes (i.e., ones obtained from Bianchi II spatially homogeneous 
spacetimes by compactifying the homogeneous spatial sections) as nilgeometric 
cosmological models (to distinguish from the conventional open models). 

By considering a compactified (spatial) manifold we have the following mer- 
its; (1) a compact Cauchy surface makes it natural to view the field equations 
as an initial value problem in cosmological context, (2) a finite spatial volume 
is physically reasonable, and (3) functional analysis on a compact manifold is 
much easier and tractable than one on an open manifold, since we can avoid 
many complexities involved in continuous spectra and those for the convergence 
of integrals (remember the contrast between Fourier series analysis on the "com- 
pactified space" = R^/Z and Fourier analysis on the open space R). 

We also mention that a motivation of considering compactified manifolds 
also comes from recent evidences that the spatial topology and global dynam- 
ical properties of solutions of Einstein's equation are related in some ways. In 
particular, recent results (l|(|5j suggest a very general picture of how an appropri- 
ately conformally transformed spatial manifold evolves in time by the vacuum 
Einstein equation, depending upon its topology. This picture motivates us to 
study linear perturbation equations for locally homogeneous solutions which 
have various spatial topologies. 

Fortunately, although it is not as trivial as the torus compactification of a 
flat manifold, the compactification of Bianchi II manifolds (or Bianchi II type 
spacetimes) is not very difficult. We will describe our compactification following 

[laiiiiTni. 

In this paper we explore detailed properties for the scalar and vector har- 
monics. We also demonstrate separation of variables for the Klein-Gordon scalar 
field equation and the source-free Maxwell equation. Although generalization 
to tensor harmonics is straightforward we leave their explicit presentation to 
a subsequent paper as well as a study of linear perturbations. This paper is 
intended to lay a solid basis for exploring those more complicated problems, 
or to be useful for many applications on the nilgeometric spacetime model like 
quantum field analysis. 

Some of the mathematical background assumed in this paper and some re- 
lated results obtained are the following. As well known (e.g., j^), the harmonics 
on a manifold AI on which a transformation group G acts are naturally obtained 
through irreducible decomposition of a representation T, called the regular rep- 
resentation, of G on L'^{M). Let g e G and let f{x) £ L^{M). The (right) 
regular representation {T,L'^{M)) is given by the homomorphism T : g ^ Tg, 
where Tg is the right translation map Tg : f{x) —> f(xa).fl^ In fact, letting 
g,g' G G, we can see TgTg,f{x) = {Tg, f){xg) = f{xgg') = Tgg>f{x), showing T 
is a representation (homomorphism), TgTgi — Tggi. This representation is how- 
ever not irreducible in general. The appropriate harmonics on M are naturally 
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obtained through an hrcducible decomposition of {T, L'^{M)). 

In our context, the transformation group G is the Bianchi II group Gn, which 
acts on M = R"^ from the left simply transitively, i.e., for arbitrary p,q G M 
there exist a unique element g G G such that gp — q. Thanks to this property, 
choosing an arbitrary fixed point o G M, e.g., the coordinate origin, one can 
identify the group Gu and the manifold M, Gu = M, by associating go with 
g. With this identification, the manifold M is also the group Gu, and therefore 
the right action of Gu on M is also naturally defined. We make M compact 
by identifying points by left action of a discrete subgroup A C Gu, M = A\M. 
The (right) regular representation (T, L'^{A\M)) on this space is indeed well de- 
fined, since keeping in mind the fact that we can identify an arbitrary function 
on A\M with an "automorphic function" f{x) on M such that f{x) = f{Ax), 
we can confirm the consistency Tgf(Ax) = f{Axg) = f{xg) ~ Tgfix). This 
shows the consistency we choose the "right" regular representation, i.e., since 
we want to define the (Killing) symmetry of the manifold with respect to the 
left action, and make a quotient by the left one, the regular representation 
on the quotient should be the right one for commutativity. The universal 
covering manifold {M,q^'^^) with a standard left invariant metric q^'^^ natu- 
rally defines a left invariant measure d/io = d^-h on M, for which we can 
define the natural inner product on L'^{M). Since for the Bianchi II group 
this measure is also right invariant d/iL oc d/XR, i.e., Gu is unimodular, the 
right regular representation (T, L'^(A\AI)) with the inner product is unitary; 
Im \fi^9)\^dno = Jm \f{x)\^dno- 

Our mode functions on M will be denoted as (/>i,m,rao (or ^i,m,no on the 
spacetime M x R) for generic modes with m ^ 0. The index m labels inequiv- 
alent irreducible representations, while the index hq labels ones in equivalent 
representations. For fixed m and no, the functions {</';};^o '^ork as a set of basis 
functions for the irreducible representation space specified by m and hq. Each 
single (f>i spans the eigenspace of an operator denoted as L^, which is like a total 
angular momentum operator. 

For the purpose of separation of variables, the most important relations 
are those for the differential representation, which is a linear transformation 
acting on the representation space spanned by mode functions. Those rela- 
tions are written in terms of group-invariant differential operators denoted as 
Xi (/ = 1 ~ 3). (See Eas. H4.14l) and (17.91) .) Indeed, since group-invariant field 
equations like the Klein-Gordon equation can be written with no explicit coor- 
dinate dependences if it is written with the invariant operators, these relations 
are found to provide the key to separate the equations. 

The formulas for the differential representation however do not provide com- 
plete information about the representation. For example, they do not tell 
whether the representation specified by m does exist in (T, L^(M)), or how 
many copies of an equivalent representation exist in (T, L^(M)). To see how 
the representation {T, L'^{M)) is decomposed to irreducible representations we 
need to find all the appropriate mode functions on M . In this paper the universal 
covering manifold M (with a group invariant standard metric) is compactified 
to a circle bundle over the torus. The regular representation (T, L^(A/)) is com- 
pletely reducible, and as a result of finding of the mode functions on the given 
topology of M we find the following: 

Theorem 1.1 Let M be the circle bundle over the torus with Euler class e — 1 
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(see Eg. 1^2.15}) for its fundamental group), and let (T, L'^{M)) be the regular 
representation of the Bianchi II group Gu. Then, it holds that 



where T„i is an infinite dimensional irreducible representation, lki,k2 a one 
dimensional irreducible representation, and the coefficient \m\ stands for the 
multiplicity in T„i. Z\{0} represents nonzero integers. 

(See 0, §10 and §11 and references therein for related mathematical works.) 
This decomposition expresses the completeness of the harmonics we construct. 
In particular, this decomposition does not depend on the Teichmiiller (or mod- 
uli) parameters of M. From the Stone-von Neumann theorem [51 122|. Tm is 
equivalent to a corresponding Schrodinger representation. 

Construction of vector (or tensor) harmonics is not difhcult on one hand. 
However the important point is to divide each irreducible space of vectors into 
subspaces such that each subspace is invariant under the action of the operator 
L^. This feature is necessary to obtain decoupled systems of ODEs when the 
background spacetime has an additional symmetry. We define three kinds of 
vector harmonics, two of which have this property. 

The plan of the paper is as follows. In the next section we describe the 
background solution and also account for some basic facts. Section |3 sets up 
some basic eigenvalues used to label mode functions, based on the compactifi- 
cation of the spatial manifold. In §^we make algebraic discussions to derive 
the %-relations. Section is devoted to construction of the mode functions on 
the spatial manifold. In section the mode functions constructed on the spa- 
tial manifold are generalized to those on the spacetime. While §§EHnideal with 
the generic modes, §0 deals with the exceptional modes, the [/(l)-symmetric 
modes, which complete all possible (scalar) modes. The results so far are applied 
to the Klein-Gordon equation and the reduced ODEs are explicitly given, with 
an asymptotic solution for a special case, in §|S1 In §|^we develop the vector 
harmonics. Section EH is devoted to an application to Maxwell's equation. The 
final section is devoted to conclusion. 

This paper is a full account, with much generalizations and development, of 
the subject outlined in §3 of 16 . Although most notations remain the same, 
one of the changes is that a quotient by the left action is now written A\AI 
instead of M /A to make clear which action is used. We employ the abstract 
index notation 24 and use leading Latin letters a, 6, • • • to denote abstract 
indices for vectors and tensors in SSlHV llOl In the other sections however we 
write them without abstract indices. We often drop the tensor product symbol 
and write, e.g., (cr^)^ instead of cr^ (g) cr^. Beware that since vectors are also 
used as differential operators, products of them like {xiY can stand for second 
order derivatives like XiXi or tensor products like Xi®Xi^ depending upon the 
quantity considered. 

2 The background solution 

Our background solution is specified by the following: (1) it is a solution of the 
vacuum Einstein equation, (2) it is spatially locally homogeneous of Bianchi II 




(1.1) 
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type, and (3) its spatial manifold is compact without boundary, in other words, 
closed. (An explicit topology will be chosen later.) 

A Bianchi type II solution is characterized by the fact that the solution (or 
the universal cover of it) is invariant under the action of the Bianchi II group 
Gii, which is a three-dimensional nilpotent (e.g., Lie group. The group 

multiplication is given by 

(a,6,c)(a',6',c') = {a + a ,b + b' , c + c + ab'), (2.1) 

for (a, b, c), (a', c') £ Gn. (Note: To save space we try to express components 
of group elements in a row form as above, but a column form is also equally 
used when it is more convenient.) 

Let M = R'^ be the simply-connected open manifold with coordinates (a:, y, z). 
We can define the group action on this manifold identifying the group manifold 
Gii with M. The left-action is therefore expressed as 

(a, b, c){x, y, z) = {a + x,b + y,c + z + ay), (2.2) 

where a = (a, b, c) € Gu, and x = {x, y, z) e M(~ Gn). Let ^/ (/ = 1, 2, 3) be 
the generators of the one-parameter subgroups (a, 0,0), (0,6,0) and (0,0, c) € 
Gii. It is easy to find they are expressed 

d d ^ d ^ d ,„ „, 

Similarly, the generators of the right actions and their dual one-forms are given 

by 

d d d d 

dx' ^'^ dy ^ dz' dz^ (2.4) 

<7^ — dx, — dy, — dz ~ xdy. 

These vectors xi and one-forms (7=1^3) are called the invariant vectors 
or one- forms of Gn, since they are left invariant; 

^«.XJ = K/,Xj]=0 = /:oa'', /,./-1^3, (2.5) 

where C^j is the Lie derivative with respect to ^/ . The invariant vectors satisfy 
the following commutation relations: 

[xi,X2]=X3, [X2,X3] = 0, [x3,Xi]=0. (2.6) 

The vectors ^/ (/ = 1,2,3) are Killing vectors for the metric of the form q ~ 
qijcr^ ® (y^ with the components g/j being constants. Riemannian manifold 
(M, q) is called homogeneous, since Gn acts transitively on it as its isometry 
group. 

A homogeneous metric is called locally rotationally symmetric (LRS) if it 
has a fourth independent Killing vector ^4. Bianchi type II LRS metrics are 
given by the metrics of the form gC^^^) = gii((cr^)^ -|- (cr^)^) -I- qssia"^)^, since 
such a metric has an additional Killing vector, given by 

d 9 1 , , 9 ,„ „, 
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This Kilhng vector generates the following one-parameter isometrics sg = e^^'* 
for the metric q^'^^^) : 




se- \ y \ ^\ ^'\y J I , (2. 

z + (g{x,y) 

1 / \ 7i"> n • 1 • . ■ r, ( COS Q — sin d \ , 

where x, y,z)£M,Kg is the rotation matrix Kg ~ \ . „ „ , and 

\^ sm fc* cos J 

Ce{x,y) = i((a;2 - y^) cos 6* - 2a;2/sin6') sin6'. (2.9) 

An LRS manifold X = (M,^'^^^)) has as a result a four dimensional isometry 
group IsomX. Let IsompX be its identity component, i.e., the component con- 
nected to the identity. An element a € IsomoA" can be uniquely expressed as 
the composite (ai, a2, 03) o sg for a choice of (ai, a2, a^) £ Gu and sg. 

The one-parameter diffeomorphism sg plays an important role even when 
the metric is not LRS. It forms a one-parameter subgroup of the automorphism 
group of Gii, which induces on the cotangent space a rotation of the invariant 
one-forms: 



sg.: I a2 j ^ | \ ) | . (2.10) 

The induced map on tangent space acts on the invariant vectors xi the same 
way with replacement Rg — *■ R-g above. The significance of the automorphisms 
of Bianchi groups was first fully recognized by Jantzen P|. Maps sg* or Sg will 
be used in this paper in several contexts. 

We can obtain the spacetime metric for the conventional Bianchi cosmology 
assuming that all the components with respect to the invariant frame (dt,(j^) 
formed by the invariant one- forms and the timelike basis dt are functions of time 
t only. The vacuum Bianchi II solution g was first obtained by Taub 1^. We 
write that metric in the following form using our invariant one-forms (|2.4|l : 

g = ^N\t)df + qi{t){a'f + q2{t){a'f + qsma^f, (2.11) 

where 

= l + l3h^P^, qi = t^P'N^, q2 = t^^^N^, 93 = IGplph^PyN^. (2.12) 

Parameters pi{i = 1, 2, 3) and f3 are constants such that P > 0, ps ^ 0, and 

Ep, = = 1. (2.13) 

When pi — P2, the solution is LRS. Although there exist two possible such 
cases (pi,P2,P3) = (0,0,1) (case I LRS) and (2/3,2/3,-1/3) (case II LRS), 
these two solutions represent equivalent one-parameter solutions. In fact, we 
can check that the case I LRS solution with P — f3i is isometric to the case II 
LRS solution with /3 = Ai = 3-2/3^-i/3_ y^^ie n we are interested in an LRS 
solution, the case II LRS solution may be preferable, since the time coordinate 
t in this solution approaches the proper time r at future infinity. This will 
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make comparisons with other models hke the Bianchi type III [201 more 
straightforward. 

It is also worth pointing out that a solution g with (pi,p2jP3) is isometric 
to another solution g' with pi and p2 swapped. In fact, it is at once using 
Ea. (|2.1()|l to see that g' is the metric induced by g' — S7r/2*g- We can 

therefore without loss of generality assume, e.g., pi < P2- 

We denote the conventional solution described so far as (M x R, g) , and 
call it the universal covering solution. On the other hand our spatially closed 
solution, denoted (M x R, g), is obtained introducing a spatial compactification 
with it. We express the solution as 

(A/xR,g) = r\(A#xR,g), (2.14) 

using an appropriate discrete subgroup T of Gn which acts spatially from the 
left on the solution. The metric g here is the one induced from g. (g and g 
are therefore locally isometric to each other.) While there are infinitely many 
possible compactifications (i.e., spatial topologies), we for definiteness specify 
the spatial manifold M to be "the circle bundle over the 2-torus with Euler 
class e = 1." (See, e.g., (3^3- In general, a closed Bianchi II manifold is a 
Seifert fiber space over a Euclidean orbifold.) The fundamental group can be 
represented in the standard notation as 

TTi{M) = (51,32,53; [51,52] = 93, [51,53] = 1, [52,53] = 1) , (2.15) 

where the brackets stand for group commutators, [a,fe] = aba^^b^^. The pro- 
cedure for the actual compactification is described in the next section. The re- 
sulting spatially compactified generalization was first constructed and discussed 
in . Note that as a result of the compactification the spatial manifold spec- 
ified by t — constant is now locally homogeneous (e.g., ^31), and the spacetime 
solution is said to be spatially locally homogeneous. 

As shown in F contains four free parameters (see Eg. HG.lll ). Our spa- 
tially closed solution H2.14|l therefore forms a six parameter solution (since the 
universal cover has as we have seen two independent parameters, f3 and, e.g.. 
Pa)- 



3 Compactification and eigenvalues 

To proceed, we need to describe the compactification of the spatial manifold 
and thereby define some eigenvalues. 

Let us first describe the canonical way of expressing a Bianchi type II locally 
homogeneous manifold (M, q) . Let g^*'-' be the standard metric given by 

= (ai)2 + (a2)2 + (a^ (3.1) 

and let 

AA= (M,e2"q(0)) (3.2) 

be the standard conformal manifold, with e^" being a constant conformal fac- 
tor. Then jOj, the manifold {M,qab) can be expressed as a quotient of such a 
standard conformal manifold (with an appropriate choice of the factor e^"), 

{M,q)=A\Af, (3.3) 
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where A is an appropriate discrete subgroup of the isomctry group of J\f, A C 
IsomA/". 

The subgroup A must be isomorphic to the fundamental group tti (M) given 
by H2.15|l . (In fact, M is a Haken manifold This means that A must 

be an embedding of tti (M) into the isometry group Isom/V. Let e IsomA/" 
(i = 1,2,3) be the image of TTi-generator gi by such an embedding. Following 
the procedure shown in miUHlj we find that it is possible to parameterize them 
in the following way: 

ai = {u, 5, 0), aa = (0, 2ttv, 0), ag = (0, 0, 2ttuv). (3.4) 

Here, a^ e Gu C IsomA/". |23 We denote A — {ai, a2, aa}. The three real param- 
eters u, V, and S are called the Teichmiiller parameters of geometric structure 
for the locally homogeneous 3-manifold {M, q) . 

We construct the harmonics on (Af, q) by two steps; first we do the construc- 
tion on a covering manifold, denoted (M, g), and then superpose appropriate 
subset of the harmonics to obtain those on {M,q). The auxiliary manifold 
{M,q) is simply defined by removing ai from A^ i.e., {M,q) = ^\A/", where 
A = {a2,a3}. ^ is a commutative subgroup of A, and therefore M is homeo- 
morphic to the much tractable manifold x R. In fact, we can easily see that 
it is each x = constant plane in Af that is compactified to a 2-torus. 

Now, let us introduce some important operators. We define the "total an- 
gular momentum-like" operator 

L^^{c^,r + {c^.,r, (3.5) 

using the Lie derivatives C^-^ and C-^^ . Note that when acting on a scalar, it 
becomes a simple form 

= (Xi)' + (X2)^ (3.6) 

which also coincides with the Laplacian Aq with respect to the standard metric 
g*-"-*, up to square of xa- 

Ao = L2 + (X3)^ (3.7) 

when acting on a scalar. Here, xi = lj2,3) are regarded as differential 
operators. It is quite important to recognize that operators xi ^-re well de- 
fined not only on the universal cover M but also on the compactified manifold 
M = A\M. In other words, the induced vector fields tt*xi on M for the covering 
map TT : M ^ M — A\M is well defined because of the invariance of xi under 
the action of yl C Gu. (For simplicity we do not explicitly write tt*, and identify 
TT*xi and xi-) Since A is also a subgroup of Gu, xi are well defined on M, also. 
The globally-defined invariant operators xi work as the fundamental derivative 
operators, since any group-invariant field equations become independent of co- 
ordinates when they are written with xi- This coordinate- free property of the 
field equations is necessary to be able to reduce the field equation to ordinary 
differential equations. 

The operator X3, called the fiber generator, has a special importance, since 
it commutes with all the invariant operators xi- In other words, xa is the center 
of the Bianchi II algebra. As a direct consequence of Schur's lemma, such an 
operator must be diagonalized to obtain an irreducible representation of the 
regular representation mentioned in Introduction. 
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The operator ~ ixi)^ + (X2)^ commutes with X3, since X3 is the center. 
We may therefore be able to diagonahze our mode functions with respect to 
both X3 and L^. Also, consider another operator ^2 = d/dy, which we can find 
commutes with both and X3, so we may diagonalize the mode functions with 
respect to ^2, also. This operator however is not well defined on M, but on 
M. This is the reason we consider the auxiliary manifold M. The existence of 
^2 is important to make it possible to perform separation of variables for the 
eigenvalue equation for (see below). 

Let ifi and be eigenvalues for the operators xs and ^2: 

X3(l> = ifJ.4>, = (3-8) 

The function cj) is supposed to be an appropriate mode function on M. Also, 
we define A by L^^ = —X^cj). The solution of these equations is given by 

4, = X{x)e"''e'''\ (3.9) 

where the function X{x) is a solution of the following "harmonic oscillator 
Schrodinger equation" : 

^ + (A^ - {fix + vf)X - 0. (3.10) 

For (j) to be well defined on M , it must be an automorphic function such that 
^(j4x) = ^(x).[2Hl From this condition we find 

fi{m) ^fn/{uv), meZ, 
v ~ i'{n) ~ njv^ n ^Tj. 

We call /i, TO, and n, respectively, the fiber eigenvalue^ fiber index, auxiliary 
eigenvalue, and auxiliary index. We call A^ the total eigenvalue. The spectrum 
of A^ is determined in the next section. 



4 Irreducible differential representations 

It is not difficult to determine the irreducible representations of the regular 
representation in their differential form, i.e., the differential representations of 
the Bianchi II algebra. In fact, it will be found that this procedure is similar to 
the one in determining quantum states of the harmonic oscillator, since (or 
the scalar Laplacian Aq) has essentially the same algebraic structure as that of 
the Hamiltonian of the oscillator. 

As mentioned in the previous section since the fiber generator xs must be 
a constant when acting on an irreducible subspace, the fiber index m € Z does 
not change values in this space. We therefore assume that to is fixed throughout 
this section. It may be helpful to bear in mind that as for the correspondence to 
quantum mechanics, the fiber eigenvalue fi = m/uv corresponds to the Planck 
constant h, while xi and X2 correspond, respectively, to the position x and 
momentum p operators, as in [x,p] = ih [xi;X2] = X3 = Remember 
however that our representation space is (M) instead of (R) . 

In this section we deal with the generic to 7^ case. The exceptional to = 
case wiU be discussed in §[7| 
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Let (j) be an eigenfunction on AI for the operators X3 L^, i.e., xs'f' — 
and L^(f> = —X^cj). It is helpful to introduce a symbol signifying the sign of the 
fiber eigenvalue, which allows us to discuss both m ^ cases simultaneously; 
we define 

<; = sign(m) = sign(/i). (m ^ 0) (4.1) 

Let us then define 

A = + '^*^2), A2 = -^iXi - <^»X2), A3 = -<iix3- (4.2) 

Then, we immediately find the following commutation relations 

[L2, A] = ^2tAiX3, [L^ A] - -^2^^2X3■ (4.3) 

This means that Ai and are, respectively, a raising and lowering operator 
for the total eigenvalue A^. In fact, since 

LMi0 = ([L2,^i] +^iL2)0 = (<^2^^lX3 + AL2)0 = -(2|/i| + A^)^^/,, (4.4) 

.4i(/) is an eigenfunction for A'^ = A^ + 2|/i|. Similarly, A2(f> is an eigenfunction 
for A'2 = A2-2|/i|. 

Taking into account the fact that Ai and A2 change the eigenvalue A^ by 
±2|/i(m)|, we can without loss of generality assume the form of spectrum as 

A^ = |/x|(2? + c™), (4.5) 

where 

; = 0, 1,---. (4.6) 

The value for / = 0, A^ = Aq = |/i|cm, corresponds to the smallest one for 
given m, which must exist because minus the Laplacian — Aq = — — (xs)^ — 
— (L^ + can have only nonnegative eigenvalues. We call I the spin index. 

At this point we know the eigenmode is specified by the pair of integers 
{l,m), so the corresponding eigenfunction can be expressed with these labels 
(pi^m- As we remarked since the value of m does not change in an irreducible 
space, we drop m and write for simplicity. 

Let us write down the whole relations we have as 

L^i = -1^1(2/ + €,„)(/),, 
A3<j)l = \n\<j)i, 
Ai(j)i = ai(j)i+i, 
-420; = Pi4>i-i, 

where we have introduced unknown constants ai and Pi , which possibly depend 
on /. 

Note the following identity that can be easily checked by a direct computa- 
tion: 

= 2^i^2 - A3. (4.8) 
Using Eas. (|4.7f) . this imphes —\fi\{2l + c„i) = 2ai-iPi — i.e., 

a,_iA=-H(Z + ^V^). (4.9) 
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Because of the arbitrariness of constant multipliers for the eigenfunctions 4>i , we 
may set ai or Pi arbitrarily, but once it is set, the other is constrained from this 
relation. We takej55j 

ai = ~-V\JA, (3i^^\{l + 2l^), (4.10) 

These do satisfy Ea. H4.9|l . Then, since we defined I so that I = gives, for given 
m ^ 0, the smallest eigenvalue of — L^, we should have 

^s'/'o = /5o0-i - v4a^(^^V^)0-i = 0, (4.11) 



implying 

So, now we have 



(4.12) 



ai^-^\fi\, A = VlAiM- (4.13) 
Gathering Eas. H4.2(I . H4.7|l . and H4.13|l . we arrive at the following set of relations: 



X2<l^i^^isj^-f{<l^i+i+l4>i-i), (4-14) 
^-0,1,- 



In particular. 



{?4>i = ~>h4>u \^i=W2l + l). (4.15) 



Now, we have found the following. A "ground state" 0o is determined as a 
solution for the two equations 

Azc^n = lAi|0o, A2c^n = 0. (4.16) 

Note that the function i^o obtained this way is automatically an eigenfunction 
of as seen from the identity (|4.8() . The excited states (t>i are determined by 
successively multiplying the raising operator (— l/-y/|/i[)yli, i.e.. 



(4.17) 



The space spanned by these functions, L^^{M) = {^"^q ai(f>i\a,i S C} n L'^{M), 
gives an irreducible subspace of L'^{M). In other words, the restriction of the 
regular representation T to L'^^{M), denoted as {Tm, L'^{M)), gives an irre- 
ducible representation. The differential representation {dTm, L'^{M)) is given 
by Eas. H4.14|l . which will be repeatedly used to separate the field equations. For 
convenience, we call these relations the x-^^^^tions. 
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5 Mode functions on the compactified manifold 



In this section we solve the eigenvahie equations for the mode functions (f)i 
under the appropriate automorphic conditions. As a result we find how many 
equivalent copies of the irreducible representation Tm are contained in T, in 
other words, the multiplicity of Tm is determined. The explicit form of itself 
is also of great interest. In this section we continue to assume m 7^ 0. 

One of the possible procedures to find explicit form of (j>i is to solve the 
equations (|4.16l) to find (/)o and compute successive differential operations in 
Eq. (|4.17|l to find general (pi . Another procedure is to directly solve the eigenvalue 
equation L^0/ — —\^(f)i for general spin index I. While both ways are possible, 
we take the latter, which provides quicker way of identifying the solutions with 
known functions. 

As remarked in §13 let us find the mode functions on M first. Note that 
Ea. (|3.10|l becomes (attaching index I to X) 

^,(,,^4)X, = 0, (5.1, 

if we define 




(5.2) 

Independent solutions to the above equation are given by Di{C,) and I?_i_i(iC), 
where Di{(^) is the Weber parabolic cylinder function. When I is zero or a 
positive integer, I?i(C) can be expressed using the Hermite polynomial i?i(C); 

Di{C)^e--^<"Hi{C). (5.3) 

Our convention for the Hermite polynomial is -ff/(C) — (— l)'e5? {S / dC})e~'^'^ . 

Since must be annihilated by A2 the appropriate choice is found to 
be Di{Q, i.e., we must take Xi — constant x Di{(). In fact, the equation 
^2^0 = together with the separation form (pi = Xi(a;)e"'^e'''^, implies {d/d( + 
(l/2)C)Xo = 0, with the solution being Xq — constant x e~*^^/'*-''» . This coin- 
cides with the one claimed for I — 0. (Conversely, as we will see, functions (pi 
constructed using these Xi can satisfy the desired relations (|4.14|l for all I, which 
justifies our choice.) 

The mode functions on M are therefore, attaching indices m and n, given 

by 




(PLm,n{^) = CiDi{±. —{fix + :^))e''^'e'''y , (5.4) 



where Ci are constants and fj, and 1/ are defined in Eas. (|3.11|l . 

The constants C/ are determined by requiring that the functions 0;,m.n(x) 
obey the x-relations (|4.14|l . Using the widely known formulas 

A'(C) = -^(A+i(C)-^A-i(C)), (5 5^ 

CA(C) = A+i(C) + ^A-i(C), 

we can easily find Ci+i — Ci, i.e., they are constants that do not depend on Z; 

Ci = C. (5.6) 
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(Actually, this is the reason we chose Eas. (|4.13|l .) The constant C may be 
determined by a normalization of the square integral on M . (See below.) 

The mode functions on M are, as mentioned, expressed as an infinite sum 
of these eigenfunctions on M . Remember that they must be invariant under 
the action oi A = {ai,a2,a3}, and the functions (j)i,m,n are already invariant 
under A = {a2,a3}. We therefore want to make a linear combination of (j)i^rn,n 
so that it is invariant under ai. Recalling the multiplication rule (|2.2fl . we find 
the following transformation law (cf. 16 , Eq.(3.15)) 

'^i,m,ri(aix) = e*""^/,„i,„+„,.(x). (5.7) 

From this we can see that the following function 4>i,m,nai defined as an infinite 
sum, is actually invariant under the action of ai (cf. 'T^, Theorem 3.1), i.e., 

oo 

<^i,™,»o(x)= e''^"°''""^^<^i,™,no+™fc(x), (5.8) 

fc— — oo 

where Z = 0, 1, • • • , oo, |m| = 1, 2, • • • , oo, no = 0, 1, • • • , |m| — 1. The sum 
is convergent at any point x. It is easy to see that since the functions 0i,m,n(x) 
satisfy the relations 14.14|l . so do (/>i.m^„Q (x). These functions are therefore the 
right mode functions on M = A\M. 

As a result of the compactification, the index no for the mode functions on 
M is now bounded by |m|. |m| is the multiplicity of the modes specified by 
the same I and m. Since for each no the functions 4>i^m span an irreducible 
subspace of L'^{M), \m\ is also the multiplicity of the irreducible representation 
Tm contained in T. 

We can summarize the results as follows. 

Theorem 5.1 There exist \m\ different sets of mode functions {'^i(^)}i^o ^^"^^ 
satisfy the relations i4-^4\ l on the compactified manifold {M^Qab) = A\Af (with 
M being the -bundle over the 2-torus with Euler class e — 1). 

Let us discuss how the mode functions can be normalized. We define the 
inner product in L^(M) as 

(/,.g)= / fg*dtio, (5.9) 

where g* is the complex conjugate of g, d^o — /\ f\ a'^ = dxdydz is the 
standard invariant measure. We want to determine the square norm 

Ni = i(bi,^i). (5.10) 

We first observe the following. 

Lemma 5.2 The invariant operators x/ (^ = 1 3j m L'^{M) are anti- 
self adjoint, x\ = -XI- 

Proof. Since 

{Xif,9)^ f {xi{f9*)-fXi9*)dfio^Ii-{f,Xi9), (5.11) 
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we need to show X/ = Jj^j xiifg*)d^J-o = 0. In fact, when, e.g., / = 1, we can 
show Ii = Jj^j Xiifg*)'^^ Aa"^ Aa'^ ^ /^^ d{fg*a'^ A a^), which is from Stoke's 
theorem Jg^^j fg*<^^ /\ = 0. Here, we have used the identity 

df - ixifW + {X2fy + (X3/)a^ (5.12) 

which is vahd for an arbitrary function / on M, and also used the relation 
(i(cr^ A cr^) = 0, which is confirmed from the definition l|2.4(l . The other cases 
/ = 2, 3 arc the same, since d{a^ A a^) = d{a^ A a^) = 0. ■ 
Remark. Operators {l/i)xi (/ = 1 ^ 3) are selfadjoint. 

Corollary 5.3 In L'^{M), A\ = -A^- 

Proof. A\ - 2-^l''{x.i + = 2-i/2(;^t _ ^,^t) ^ _2-i/2(^, _ ^^^^) ^ 

-A2. ■ 

Returning to the issue of A^j, consider A^;+i = When m 7^ 0, 

from Eas. l|4.7|l and Corollary 15. 31 we have 

iVi+l = ^--^(^i</);,^i</);) 

= ^^(0i,-42^i</'O (5.13) 
a; a; 



-A+i 



OLl 

Substituting our choice H4.13|l of a; and (3i we have 

Ni+i = {l + l)Ni. (5.14) 

Taking A'o = 1 we conclude 

Ni = II. (5.15) 

Now, we have the following. 

Theorem 5.4 Suppose that (j)i are mode functions on {M, qat) — A\J\f such 
that they satisfy the relations j^. J^l ^- Multiplying the same constant normaliza- 
tion factor C to all 4)1, (pi — > C(/)i, does not change those relations. By choosing 
C appropriately, we can make the normalization 

{4'Lm,no,4>l'.Tn',n'J = l^- Sw Smm' Sngn'^- (5.16) 

hold. 

Proof. The orthogonality for m and m' is apparent from the fact that 
/i(m) is the eigenvalue of the selfadjoint operator {l/i)x3- The orthogonality 
for hq and u'q comes from the orthogonality among the mode functions on M: 

{4>l,m,n,4>l,m,n')M = _ 4'l,m.,n(t>lm,n'dtiO = 0, {nj^n'), (5-17) 
JM 
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which is also apparent from the fact that the operator (1/1)^2, of which eigenval- 
ues are vin), is selfadjoint on L^(M). Observing that when uq ^ n^, 4'i,m,nQ and 
4'i,m,n'„ are linear combinations in different sets {4>i,m,n}neAri and {4>i,m,n}neAf2y 
A/i n A/2 = 0, we can easily see 

/ (f>l,m,no(l)l„,^n'df^O = 0, ("O 7^ "o) (5.18) 

Jm 

which in tern implies the orthogonality for uq and Uq in L'^{M). The other part 
has already been proven. ■ 
Remark. The constant C does not depend on /, but can depend on m and 
no, so we may write C = Cm,no- 

Apparently, if we define 

<^f"Uo = (5.19) 
they become orthonormal to each other: 

(4m.«o' 4"L,n;,) = 5ll' 5mm' 5n„n;- (5.20) 

As seen from Eas. (|5.13|l and H4.9|l this corresponds to choosing 

a|") =e^«Vl/^l(^ + l), /^|"^=-e-*®'^, (5.21) 

where 0/ is an arbitrary phase factor, which we may want to take zero, Qi = 0. 
Substituting ai = a|"' and /?; = into Eas. (|4.7|l we have another version of 

X" relations for <j)f]l^ „^ , which have the most direct correspondence to the usual 
relations between the quantum states of the harmonic oscillator. As mentioned 
however we employ the unnormalized (pi in this paper for the convenience of 
computations. 



6 Further transformation 

The mode functions 0i,m,no(x) shown in the previous section are not well defined 
on the general spacetime solution [M x R, g) = T\{M x R, g), because of the 
fact A ^T. Remember that T C Gu is a four-parameter embedding, while A is 
a three-parameter one. This incompatibility means that we cannot identify the 
coordinates x in A/" with the spatial coordinates x in [M x R, g) . [SIl Although 
it is expected that an appropriate diffeomorphism can make the mode functions 
well defined on the spacetime, such a diffeomorphism can affect the x-relations. 
In the following, we show by explicit computations that this is the case but a 
further renormalization makes the mode functions retain the original x-relations. 
As shown in |18j , the covering group T can be parameterized as 



r = {71,72,73} = 




15 



where 7i E Gu, and 9, u, 5, and v are real parameters. 12] EH Let {M,q) = 
T\{M,qab) be the spatial section of the solution, where q is the spatial part of 
g. For convenience of considering diffeomorphisms between the spatial universal 
cover {M,qab) and the standard conformal manifold A/" = (M, e^"q^°^), let us 
distinguish the latter manifold (without metric structure) by denoting M'. This 
distinction is meant to be helpful just to keep track of the direction of the 
diffeomorphisms we consider. 
Let ip he a diffeomorphism: 

ip-.M^M' (6.2) 

such that 

r = ^p-^ o Ao^. (6.3) 

Then, the induced function of (pi^m.noi^) will be invariant under F, since so is 
(t)i,m,no (x) under A. 

We can find '0 in Isonx/V, which is explicitly given by 

V' = bose, (6.4) 

where b e Gu is 

b = {ttv cos 6* sin 6*, u^^{~(^^g{u, 5) + ttSv sin 6 cos 6),0). (6.5) 

We have, for simplicity, chosen the third component of b as zero, though it 
can be an arbitrary constant. The induced vectors ip^'xi on M' becomes a 
linear transformation of Xh due to the property that elements in IsomA/" are 
automorphisms of Gu. Note that when acting on the vectors x/j the induced 
map -0* = b* o Sg becomes the same as Sg, since xi are by definition invariant 
under the induced map b* for b G Gu- Therefore from (the vector version of) 
Ea. (|2.1U|l . we have 







/ rxi 


r : 


i)- 


- rx2 






V rx3 



R-g 

^3 



Xi 

X2 ; I . (6.6) 



Let (pl (x) = ('0*(/);)(x) = 0(O0(x), where x e M. (Superscript (ss) stands for 
"spatial section" .) Then, 

xifpf""^ = xi(V'*0i) = (V'*xi)'/'i ° 

= (cos^xi + sin 0X2)0; ° V' 



2 (-e-'Vi+i+e^'/^i-OoV' 

where we have used the relations (I4.14|) . Similarly, we obtain 

X20^^ = ^^^J^ {e-'Uf:\ + e^'ld^fll) , (6.8) 

and X30/^''^ = i^(j3f^\ These relations are different from the original x-relations 
(|4.14|) unless is a multiple of 27r. 
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However, it is possible to renormalize 0; so that the original x-i'elations 
are recovered. It is straightforward to check that 

= e^^Vr^^ (6.9) 

gives such a renormalized function. Thus, we have found that the functions 
(fi = €^^^4)1 o ip with (pi given by Ea. H5.8|l and ip being Ea. (|6.4|l . are served as 
the right mode functions on the spatial section (and therefore on the spacetime) 
that satisfy the relations H4.14|l . This provides a direct proof of the following: 

Theorem 6.1 There exist \m\ different sets of time-independent mode func- 
tions {¥'i,m(x)}^Q on the spatially closed Bianchi II solution {M x R, g) — 
r\(M X R, g) such that they satisfy the relations \4-.14^ . 



7 [/(l)-symmetric modes 



Let us, for completeness, consider the modes with the fiber index m being zero. 
We call these modes U{1)- symmetric, since they are constant along the U{1) 
(~ S^) fibers. 

Let (j) be an eigenfunction for m = 0, i.e., xs'P = {dtp/dz) = 0. This in turn 
implies that xi X2 are commutative when acting on 

[xi,X2]</' = 0, (7.1) 

since [xi:X2] = X3- -Due to this property, the harmonics describing the m — Q 
subspace of L'^{M) become the usual Fourier expansion on a torus. (Thus, 
together with the results for the generic modes we obtain Theorem II. Ih In 
the following we explicitly determine the spectrum of the eigenvalues for the 
operators xi and X2 in terms of the spacetime moduli parameters u, 6, v, and 

e. 

Let us first work on the 3-manifold A\J\f (not on the spacetime manifold) 
as we did for the generic (m ^ 0) case. Taking the form of A into account, we 
label the mode functions with the the following equations 

("Xi +(5X2)0 = 27rifci(/), ^^^^ 

where the eigenvalues ki and k2 are to be used as labels. Let us therefore write 
the solution of these equations as (p — 4if \ (x) — {x,y), which is given 

by 

0^°)^^ ^ constant x ei(2'^'^-i^+^(-'5^+"a))_ (7.3) 
Since, recalling the rule H2.2|l . 



(aix) = ix + u,y + S) = 4°),^ (x)^^-'^^ 

<fc.(a2x) = + 27r^;) = <^,Jx)e2"'=% 

we find 



(7.4) 



ki e z, fcs e z, (7.5) 



17 



for to be well defined on M ^ A\M. The remaining condition (p'f^i^^ (^ax) = 

4>k^M trivial. 

Therefore from Eas. l|7.2|l the x-relations for the C/(l)-symmetric modes are 
X2<,,-*i^2(fc2)<,,, (7.6) 



where 



Ki{ki,k2) = -{2TTki - -fca), 

I ^ (7.7) 

i^2(fc2) ^ ^. 

V 

To extend the mode functions on the spacetime manifold r\ (ill x R, g) , we 
need to apply the diffeomorphism ip defined in Ea. Ht).4|) again, and as a result 
the spectrum of the eigenvalues are altered. (Contrary to the generic case, any 
renormalizations of the resulting mode functions do not affect the x-relations.) 

As in the previous section, let us define ^l^^ = ° Then, from 

Ea. (|6.()(l . we have, e.g., 

= (cos(?xi+sin(?X2)</>Sfc,oV' (7.8) 

= i{coseKi+fiineK2)Vk!M- 

A similar result is also obtained for X^vfi k^' write the final form of the 
relations as follows. 

X2fklk2 ^ «'«2(fci, A:2)¥'^°^_fc^, (7.9) 



where 



Ki(fci, fc2) = cos6'i^i(fci, A;2) + sin6'i^2(fc2), 
K2(fci, fc2) = - si\i9Ki{ki,k2) + cos0K2{k2). 



(7.10) 



Now, we have: 

Theorem 7.1 There exist time-independent mode functions (^^^^^^(x), fci,/c2 G 

Z, on the spatially closed Bianchi II solution {M x R, g) = T\{M x R, g) such 
that they satisfy the relations \7.9^ . 

8 Application to the Klein-Gordon equation 

As an example, let us consider the Klein-Gordon equation 

(g'^^VaVb - m|)$ = 0, (8.1) 
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where ni$ > is the mass of the field $. Va is the covariant derivative operator 
associated with the spacetime metric gob- It is straightforward to see that this 
equation on our background can be expressed, using the invariant operators x/, 
as 

^ ^ ^^N-^i) + A, - ml) $ = 0, (8.2) 



where Ag is the Laplacian with respect to the spatial metric Qat', 

A, - gr'(xi)' + 92\x2r + q^\xz)\ (8.3) 



and = V-detgafc = 4|p3|/3t7V2. 

Let us consider a generic irreducible component of i.e., $ = $m,noi ™ 7^ 0. 
We can expand this component as 

oc 

$(t,x) =^a,(t)(^,(x), (8.4) 

where (/p; = (fii^m.na are the mode functions mentioned in Theorem 16. II 
From the relations 14.14|l . we have 



2,„ _ ,,2, 



(8.5) 



from which we immediately obtain the following wave equations for ai (t) : 

di + jai + Z{t)ai = J(i; a;_2, ai+2), (8.6) 



where 

16(p3)2/?2 



(8.7) 



with the inhomogeneous term / being 

/(t; ai.2,ai+2) = ^(^-2^1 _ i-2p.)^_^ ^ ^ 2)(Z + 1)0^+2). (8.8) 

(In I{t; a;_2, a(+2), ai-2 should be regarded zero when I — and 1.) 

Note that the inhomogeneous term / introduces couplings with the next 
neighboring modes with / ± 2. The equations H8.6|l therefore comprise two sys- 
tems of infinite number of equations, the one with / = even and the one with 
I — odd, unless the background is LRS. When on the other hand the back- 
ground is LRS, each equation H8.6|l for a given / becomes closed itself, due to 
the vanishing of the inhomogeneous term /. 

When the background is LRS, we can find future (i — )■ 00) asymptotic solu- 
tions: 
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Proposition 8.1 On the LRS Bianchi II vacuum solution with pi = P2 = 
2/3 and = —1/3, the scalar field equation ^8. 6)) for a generic mode has the 
following fundamental solutions as t ^ oo: 

y^{t) = t-§e±vTKG(t)(i + o(l)), (8.9) 

where 

The symbol o(l) stands for a function such that hmj^oo o(l) ~ 0. 

Proof. This result is a generalization of Theorem 3.4, Ref.^Q, with finite 
mass m$. As emphasized there, it is an appropriate choice of new time vari- 
able that is essential to obtain an asymptotic solution. In the present case an 
appropriate choice T{t) is given by 

dT 3 2/3m| /3/3 8^\ 1 

Following the procedure shown in the reference, one obtains the asymptotic 
solution H8.9(l . ■ 
It is worth noticing that the asymptotic solution only depends on the fiber 
index m and the other index I does not affect them. See the final section for 
more discussion. 



9 Vector harmonics 

Let us discuss how we can construct the vector harmonics. 

Fortunately, this is on one hand a trivial issue as the invariant frame {cr^, xi} 
is well defined on the compactified manifold M (and on the spacetime manifold 
M X R). This means we can define the components Tj...'^ " of any sort of tensor 
Ta---^ " with respect to this invariant frame; 

TiJ- ^Ta...''-xi''---(^\-- - ■ (9.1) 

All we have to do is to expand these components with respect to the scalar 
harmonics 0;,m,no and <j)f^ j.^- This procedure corresponds to using the set of 
one-forms {(pic^a, 4>i<^'^ a, 4>i'^'^ alf^Q as vector harmonics for given m ^ and uq. 
(As for m = 0, should of course be replaced by 4>^^^ k^i-) spacetime 
manifold, we also need the timelike mode vectors, as well as the replacement 
4>i,m,no ~^ ^i,m,no- We therefore can define the following harmonics (basis mode 
vectors) 

(l^'°)a = ipi{dt)a, iV{% = Vict\, {V{^)a = cpia\, (y/3)a = cpia\, (9.2) 

(or the one with y^^^^^^^ instead of ipi for m — 0). We call them the simple vector 
harmonics. 

The advantage of considering these harmonics is that they are apparently 
complete (since the scalar harmonics used are complete). However, this choice of 
harmonics does not reflect very well the group structure regarding the rotational 
automorphisms (|2.1(J|) . A more natural and convenient choice can be obtained 
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by using the "spherical bases" m] Aj (/ = 1 ~ 3) defined in Eas. (|4.2f) . or their 
duals such that Ai"'g'^a = ^j', 

= -^{a^ - ';ia^), = -^(ct^ + <iia^), g^ = <,ia^ . (9.3) 

Using these we can set up new harmonics as 

(y,°), = ^i{dt)a, {Vi^)a = Vl+ig\, (T^')a = Vl-ig\, iV{')a = ifiW^. (9.4) 

(When m — we just use the same harmonics as the simple harmonics.) Be- 
ware that (pi±i are used to define (VJ^)a and (VJ^)a. The reason will become 
clear below. Because of this index correspondence, we should think that the 
harmonics {(VJ^)a}/=o defined for I > —1 (not for I > 0). The basis (Vi;^)a 
for / = — 1 is nonzero only for I — 1, and the others should simply be regarded 
as zero. We call these harmonics the polarized vector harmonics or the standard 
vector harmonics. 

We can confirm that for given m, the simple harmonics and the polarized 
harmonics span the same space of vector fields. This ensures the completeness 
of the polarized vector harmonics. 

Theorem 9.1 For given fiber index m e Z (and given auxiliary index uq), the 
linear span of the simple harmonics 

( oo 3 



Span(F"")^0 }_^cz,(^'Oa 



CIA e C j , (to ^ 0) (9.5) 
(the m — case is defined similarly) and that of the polarized harmonics 



Span(F")^ EE^^''(^')- 



ci,i e C L (to ^ 0) (9.6) 



(the m — case is defined similarly) are the same; 

Span(T/'™) ^ Span(F™). (9.7) 

Proof. This is trivial for the m = case, since in this case the two kinds of 
harmonics are the same. So, we can assume to 7^ 0. From the definition H9.4|l . 
it is apparent that each basis one-form VJ^ of the polarized harmonics for given 
TO can be expressed as a linear combination of the simple harmonics belonging 
to the same to. Conversely, each basis one-form V/^ of the simple harmonics for 
given m can be expressed as a linear combination of the polarized harmonics 
with the same to as 



(9.8) 





= vr 
















vr 


= vr 
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Therefore the two sets are related by a regular linear transformation, which 
proves the identity of the two spans. ■ 
The significance of the polarized harmonics is that for given / (and as usual, 
given m and no), the space spanned by them 



cieC\ (9.9) 



1 7=0 ) 

is invariant under the operation of L^. To show this, let us start with observing 
the commutation relations 

[AuA2] = ^3, [AuAs] = 0, [^2,^3] = 0, (9.10) 

from which one can immediately have 

LajAj - [Ai,Aj] - e/j3-43, (9.11) 

where cjjk is the unit skew symmetric symbol; €123 = +1, eijK = ^[ijk\- Then, 
noting the duality Aj°'q^ a — 5^ : it is also easy to see 

Ca,q' ^-5ieiKzQ''. (9.12) 
From this equation and the ^-relations, as well as the identity 

= £^j£^2 + £^2'C^u (9.13) 

one obtains 

L^i^w') - -X^ivw' - 2^\5i{^i+ig' + lipi^ig"). (9.14) 

This equation happens to be valid when / = 0, as well. Converting to our bases 
(|9.4(l . we obtain the following statement. 

Lemma 9.2 Let (I — ^ 3) be the basis mode one-forms defined in 
Eqs . {9.4}^ , and let be the second order Lie derivative operator defined in 
Eq. \3. Then, it holds that 



(9.15) 



where ~Xf is the eigenvalue of with respect to the mode function ipi, defined 
in gg. 177751 ). 

Note that the right hand sides of Eas. H9.15|l are linear combinations of {VJ^}/=o 
belonging to given I. This proves the invariance we claimed: 

Theorem 9.3 The linear span Span(V;), defined in Eg. \y.y\) . of the harmonics 
{^'^}/=o invariant under the operation of : 

L2Span(V^) C Span(Vi). (9.16) 
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Thanks to this property, if the background spacetime is LRS and the field 
equation is invariant under the LRS-action induced by se as well as the usual 
group action, the ODEs reduced fi'om the field equation become all independent 
from the others. In other words, "each Z" decouples from the others in case of 
LRS. (This property does not occur for the simple harmonics, since the linear 
span of them for a given I is not L^-invariant.) 

Before ending this section let us mention the approach taken in JS], Sec. 3. 4, 
which was based on the analogy of the spherically symmetric case |13) or the 
Bianchi III hyperbolically symmetric case [2111 . Let us denote the harmonics 
used in the reference with two dashes like (V/'^)a. V/'^ and V/'^ are defined in 
the same way (up to constant multiplication factor) as the polarized ones; 

Vr^V^', Vr^,^^V^^ (9.17) 

V/'^ is called the time-like basis one-form, while V/'^ is called the fiber basis 
one- form. Let us then consider the "plane field" spanned by xi and X2, which 
is horizontal to the base manifold. This horizontal plane field is invariant under 
the natural actions of the fibers generated by X3- We define the "area two- form" 
e of this field by Sab — "^cr^ [a'^'^ b] (although the plane field is not integrable). We 
define the horizontal gradient (HG) basis one- form (corresponding to the "even" 
one-form ^11201) by taking gradient of the scalar harmonics and subtracting 
the fiber part; 



{vr)a=davi-{vr)a 



(9.18) 
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And, we define the dual horizontal gradient (DHG) basis one- form (correspond- 
ing to the "odd" one- form) by taking the dual gradient associated with e; 



= ■i{{X2(Pl)cr^ a - (xl<^/)o-^a) 



(9.19) 
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{ipi+i + lipi-l) a^a + i - IVl-l) O'^a) 



To raise an index for Sab we use the (inverse of the) standard Bianchi II met- 
ric ft.C)"'' = xi°'Xi' + X2"'X.2^ + Xs^Xs^- It is clear that the harmonics VJ" = 
{(VJ"^)a}7=o ^''^ equivalent (i.e., their span is the same) to the polarized har- 
monics Vi = {{V/)a}^^Q for each I, since they are merely related to each other 
by a regular linear transformation, as seen from the relations 

vr^-^/\fAiVl'-lV,'), vr^-^V\lA{Vi' + lvn, (9.20) 

as well as (19.17(1 . We call these harmonics the mixed vector harmonics. This 
choice is particularly convenient for Maxwell's equation, since the J7(l)-gauge 
transformation, a shifting of the vector potential by a gradient of scalar, does not 
affect the component for the DHG basis one-form. (As a result, this component 
itself becomes gauge- invariant. See the next section.) 
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10 Application to Maxwell's equation 

As an application let us consider the source-free Maxwell equation VFab = 0. 
Since the electromagnetic field Fat is given by (twice) the exterior derivative of 
the vector potential Aa', Fat — daAb — dbAa, this equation can be dealt with by 
vector harmonics. 

Let us consider the irreducible component belonging to given m ^ and no- 
For this component we can expand the vector potential as follows: 

oo 3 

^'^^ EE7fw(^/")- (10-1) 

7=0 

where (V^"^)a are the mixed vector harmonics. The four kinds of functions of 
time 7^'' (t) (/ = ~ 3) serve as the field variables. 

The quantities we are interested in are the C/(l)-gauge invariant variables, 
which can be easily found by inspecting components of the field strength Fat — 
daAb — dbAa- We obtain the following four independent f7(l)-invariant variables: 

Qi'^ - li'' - Q'i' - 7f , = i? - P^' - i^'- (10.2) 

Although function P3 = 73 — 70 is also invariant, it is found that it can be 
(consistently) solved with the others, due to the constraint part of Maxwell's 
equation = {dt)°''^^Fab, which can be written, using the invariant operator, 
as 

= N-\q^\^A^ + q^\2A2 + 93-^X3^3) - AgAo, (10.3) 

where Ao = N^^{dtY Aa, a-^^^ = Xi"'^a (/ = 1 ^ 3). Laplacian A, is given 
in Ea. (|8.3(l . The evolution equations = xi"'^^Fab (/ = 1 ^ 3) can similarly 
be written 

= -N-^F^i + q^\iFu + q^\2F2i + Qs'xsFsi + N-\Y^qiFoi 

- {2N)-\qY\, + qY'q2 + qY'q3)Foi + (<Zig2)"'(73'5?Fi2, (10.4) 
(no sum for repeated indices) 

where Fqi = N-\dtTxi''Fab and Fjj = Xi^x/Pab- After a rather lengthy 
computation they become the following: 















— ^2 ' 
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93 
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^2 




1 93 

2<Z3 




93 



(10.5) 
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where the inhomogeneous terms are: 



2|mI 



^iQi' - 92"')lMl{a + 2) (g^'^ + ^{21 + 5)Q^t' 



4 V9i 92 



+ 2) (Pi''+^' - ^Pa'' 



(1+2) 
2 



+r^(pr^)+.pr^))} 

(gr' - 92"')lMi|a + 2) (^g^" + (2^ + 5)0' 

+ ri(,gr^) + (2?-3)g^^)) 
1^91 92M^^^2)(,p('+^)-p('+2') 
-ri(,p('-2) + p('-2)' 



(10.6) 



4 \qi 92 



Again, we can see the same quahtative features in the Klein-Gordon equa- 
tions; we have two systems of infinite number of equations, the one with I = even 
and the one with I = odd, unless the background is LRS. When the background 
is LRS, the couplings between mode I and the next neighbors / ± 2 are cut off 
due to the vanishing of the inhomogeneous terms /g^ , /p^ and Ip^ and this 
makes each system of four first-order equations (|10.5() for given / closed itself. 
This is of course a result of the invariance described in Theorem IH. 31 

11 Conclusion 

There are three main results about the basic properties for the nilgeometric 
harmonics obtained in this paper. They are (i) the irreducible decomposition 
of the regular representation fTheorem ll.l|l . (ii) the explicit form of the mode 
functions, and (iii) the differential representation formula, the ^-relations (see 
Eas. H4.14|) and {^j)). The decomposition (i) represents the completeness of our 
harmonics. As for the point (ii), remember that we have two kinds of formula, 
the one for the canonical manifold and the one for the spacetime manifold. The 
former is given in Eas. H5.8| l and H7.3|l . while the latter is obtained by the trans- 
formations explicitly given in § El and §0 Remember also that the x-relations 
are the most important for the purpose of separation of variables. We also have 
generalized the (scalar) harmonics to vector harmonics, and demonstrated sep- 
aration of variables for a scalar equation (the Klein-Gordon equation) and a 
vector equation (the Maxwell equation). 
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As we have seen, when the fiber index m is nonzero the ODEs reduced 
from a field equation, e.g., the ones from the KG equation, become systems 
of infinite number of simultaneous equations. In this sense, infinite number 
of different modes are coupled to each other. This is a result of the fact that 
the corresponding irreducible representation is infinite dimensional. When the 
background is LRS however, the couplings between the modes are cut off and 
as a result, each single reduced KG equation becomes closed itself. Although 
the Maxwell equations give rise to much more complicated reduced equations 
because of the multiple components of the field variable, it has the same feature 
that the couplings between modes disappear when the background is LRS. It is 
also apparent that the linear perturbation equations will have the same feature 
if we choose the tensor harmonics so as to possess the invariance under like 
the vector harmonics do (cf. Theorem 19.3(1 . 

An interesting fact is that as shown in §|H1 the future asymptotic solution 
of the LRS KG wave equation depends only on the fiber index m and does not 
depend on the spin index I. This fact seems to indicate a clue for analyzing the 
generic non-LRS cases, since it suggests that the couplings between the modes 
asymptotically disappear even when the background is non-LRS, at least if the 
background is close enough to the LRS one. See ^7] for the same ("fiber term 
dominated" ) behavior of other models (Bianchi VIII and III) . Detailed studies 
of the non-LRS cases, as well as linear perturbations of the nilgeometric model, 
will be reported elsewhere, on the basis of this work. 
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